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ABSTRACT 

Q_i| The growth of the angular momentum L of protogalaxies induced by tidal torques is 

■ reconsidered. We adopt White's formalism and study the evolution of L in Lagrangian 

coordinates; the motion of the fluid elements is described by the Zel'dovich approxima- 
tion. We obtain a general expression for the ensemble expectation value of the square 
of L in terms of the first and second invariant of the inertia tensor of the Lagrangian 
volume r enclosing the protoobject's collapsing mass. We then specialize the formal- 
' ism to the particular case in which T is centred on a peak of the smoothed Gaussian 

density field and approximated by an isodensity ellipsoid. The result is the appropriate 
analytical estimate for the rms angular momentum of peaks to be compared against 
simulations that make use of the Hoffman-Ribak algorithm to set up a constrained 
f^S density field that contains a peak with given shape. Extending the work of Heavens 

\D . & Peacock, we calculate the joint probability distribution function for several spin 

parameters and peak mass M using the distribution of peak shapes, for different ini- 
. tial power spectra. The probability distribution for the rms final angular momentum 

Q-T (L 2 ) 1 / 2 on the scales corresponding to common bright galaxies, M « 10 11 M Q , is cen- 

Q ■ tred on a value of ss 10 kg m 2 s~ 1 , for any cosmologically relevant power spectrum, in 

^—j | line with previous theoretical and observational estimates for L f . Other astrophysical 

consequences are discussed. In particular, we find that typical values (A 2 ) 1 / 2 « 0.1 of 
the dimensionless spin parameter for peaks smoothed on galactic scales and of height 
v ~ 1, usually associated with late type galaxies, may be recovered in the framework 
of the Gaussian peak formalism. This partially relaxes the importance attributed to 
$2] ■ dissipative processes in generating such high values of centrifugal support for spiral 

galaxies. In addition, the values of the specific angular momentum versus mass - 
as deduced from observations of rotational velocities and photometric radii of spiral 
galaxies - are well fitted by our theoretical isoprobability contours. In contrast, the 
observed lower values for the specific angular momentum for ellipticals of the same 
mass cannot be accounted for within our linear regime investigation, highlighting the 
importance of strongly non-linear phenomena to explain the spin of such objects. 

Key words: galaxies: formation - large-scale structure of the Universe 



1 INTRODUCTION 

It has been argued that tidal coupling between the inhomogeneities in the primordial matter distribution, in the context 
of a gravitational hierarchical scenario, may explain the acquisition of the angular momentum by a protogalaxy. This idea, 
originally due to Hoyle (1949) and applied by Sciama (1955), has been first thoroughly examined by Peebles (1969), who 
demonstrated that the tidal spin growth of the matter contained in a spherical (Eulerian) volume is proportional to t 5//3 
in an Einstein-de Sitter universe (t is the standard cosmic time). Specifically, Peebles' analysis is based on a second-order 
perturbative description, since a spherical volume does not gain angular momentum from tidal torques in linear approximation, 
as pointed out by Zel'dovich and reported by Doroshkevich (1970). An important result of Doroshkevich's paper is that a 
generic nonspherical volume enclosing the protogalaxy acquires tidal angular momentum proportionally to the cosmic time 
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t during the linear regime. This theoretical prediction has been confirmed by the iV-body simulations of White (1984). In 
addition, White showed that the second-order growth described by Peebles is due to convective motion of matter across the 
surface of the initial volume T containing the protoobject. An important point of Peebles and White's theoretical analyses 
is that they describe the tidal torques acting on a volume centred on a random point in the smoothed density field, which 
does not necessarily enclose a bound protogalaxy. In contrast, galaxies are expected to form around (high) peaks on relevant 
scales in the density field. This idea, in embryo in Doroshkevich (1970), developed into the biased galaxy formation scenario, 
where only density maxima above a given threshold (peaks) of the initial Gaussian density field can eventually form galaxies 
(Kaiser 1984; Politzer & Wise 1984; Peacock & Heavens 1985; Bardeen et al. 1986). 

The acquisition of angular momentum due to tidal torques by Gaussian high-density peaks has been recently analyzed, 
amongst others, by Hoffman (1986; 1988), Heavens & Peacock (1988), Ryden (1988), Quinn & Binney (1992) and Eisenstein & 
Loeb (1995). Comparisons with iV-body simulations are displayed in Efstathiou & Jones (1979), Barnes & Efstathiou (1987) 
and Warren et al. (1992). 

In this paper, we re-examine the growth during the linear regime of the angular momentum L of protogalaxies induced 
by tidal couplings with the surrounding matter inhomogeneities. The layout of this article is the following. In the next 
section we briefly review White's formalism describing the linear evolution of tidal galactic spin. The motion of the mass 
fluid elements is described by the Zel'dovich approximation: the invariance of the angular momentum with respect to the 
Eulerian and Lagrangian description is stressed. Next, we derive a general (but approximate) expression for the ensemble 
expectation value of the square of L, (L 2 ), in terms of the first and second invariant of the inertia tensor of the Lagrangian 
volume r. We then specialize our formalism to the particular case in which the Lagrangian volume is centred on a peak of the 
underlying smoothed Gaussian density field and approximated by an isodensity profile ellipsoid: in this case, we obtain the 
correct constrained ensemble average for these objects with preselected inertia tensor. The result is the appropriate analytical 
estimate for the rms angular momentum of peaks to be compared against simulations that make use of the Hoffman-Ribak 
algorithm to set up a constrained density field that contains a peak with given shape (Hoffman & Ribak 1991; van de Weygaert 
& Bertschinger 1996) . Extending the work of Heavens & Peacock, we calculate the joint probability distribution function for 
several spin parameters (angular momentum, specific angular momentum, angular momentum in units of M 5 / 3 ) and peak 
mass M using the distribution of peak shapes (Bardeen et al. 1986), for different density power spectra. Finally, we discuss 
astrophysical implications in the last section. Technical details are given in appendices. 



2 ANGULAR MOMENTUM 

We carry out our analysis of the evolution of the angular momentum in three steps: first, we review White's method for 
obtaining an expression for L involving the shape of the protoobject and the distribution of the surrounding matter. Next, 
we simplify the expression for L by performing the ensemble average in subsection 2.2. Finally, we specialize the formalism 
to the case where the object is centered on a peak of the (Gaussian) underlying density distribution. 

2.1 Dynamical description 

Let us assume that at the epoch of structure formation the matter may be described on the relevant scales as a Newtonian 
collisionless cold fluid (dust) embedded in an expanding Friedmann universe with arbitrary density parameter fi (but, for 
simplicity, vanishing cosmological constant). We indicate the comoving Eulerian spatial coordinate by x. The physical distance 
is r = a(t)x, with a(t) the expansion scale factor. In the Einstein-de Sitter model, a(t) tx t 2 ^ 3 . 

The angular momentum L of the matter contained at time t in a volume V of the Eulerian x-space is 



Here, p = pb(l + S) denotes the matter density field, pb the background mean density, S the density fluctuation field; v = dr/dt 
indicates the velocity field and u = adx/dt is the so called "peculiar" velocity (see, e.g., Peebles 1980). The origin of the 
Cartesian coordinate system is assumed to coincide with the centre of mass. Since we are interested in the intrinsic angular 
motion, we disregard the centre-of-mass motion. 

It is important to note that the previous integral over the Eulerian volume V may equally be written as an integral over 
the corresponding Lagrangian volume V: 




(1) 




(2) 



where p&a 3 = po«o — Vo in the matter dominated era. Equation 



can be obtained from equation (hi) using the mapping 



x(q,i) = q+S(q,i) 



(3) 



© 0000 RAS, MNRAS 000, 000-000 



Linear evolution of the tidal angular momentum 3 



from Lagrangian coordinates q to Eulerian coordinates x, where S is the displacement vector. The determinant J of the 
Jacobian of the mapping q — > x(q, t), dx J -1 = dq, is related to the density fluctuation through the continuity equation: 

l + *[x(q, *),*] = Ji^ty 1 . (4) 

Before the occurrence of shell crossing (caustic formation process) J is non vanishing (see, e.g., Shandarin & Zel'dovich 1989). 

The expression for L(t) in equation (^) allows us to apply the Lagrangian theory directly: we stress that it is an 
exact relation. Consequently, when applying perturbation theory to equation (^), perturbative corrections to S (Bouchet et 
al. 1992; Buchert 1994; Catelan 1995 and references therein) provide perturbative corrections to L (Catelan & Theuns 1996). 
Furthermore, equation ^ shows that L(t) depends on the shape of the Lagrangian volume T, which encloses, by definition, all 
the fluid elements that eventually form the galaxy. The boundary of a given condensation region F depends on the particular 
realisation of the density field and in general it may be extremely fuzzy. 

The linear Lagrangian theory corresponds to the Zel'dovich approximation, 

S(q,t)^S (1) (q,t)=D(t)VVi(q), (5) 

which coincides with the linear Eulerian approximation to first order in the spatial positions, and gives a qualitative correct 
quasi-nonlinear solution by extrapolating the straight trajectories of equation (|^) beyond the range of strict validity (Zel'dovich 
1970a, b; Shandarin & Zel'dovich 1989). The temporal function D(t) describes the growing mode of the density perturbations: 
in the Einstein-de Sitter model D(t) oc a(t). We neglect the decaying mode. 

The potential ipi(q) is the linear gravitational potential (see the discussion in Kofman 1991): if i/>i(p) = J dqV>i(q) e~ lp q 
is the Fourier transform of the gravitational potential, then the Fourier transform <5i(p) of the linear density field 5^'(q, t) = 
D(t)Si(q) is related to V'i(p) via the Poisson equation, <5i(p) = p 2 ?/>i(p), where p = |p|, and p is the comoving Lagrangian 
wavevector. 

In the Zel'dovich approximation, the exact equation (0) simplifies to: 



(6) 



L(t) » L (1) (t) =»a 2 | dq (q + S (1) ) x = Vo a 2 ^ ^ dqq x V*(q) . 

As in White (1984), if it is assumed that -0i (q) is adequately represented in the generic volume F by the first three terms 
of the Taylor series about the origin q = 0, then each component L a l \t) may be written in a compact form as (White 1984) 

Z#> (i) = a(t) 2 D(t) e a ^ 23« , (7) 

once the deformation tensor at the origin q = 

V$=V^(O) = d d„MO), (8) 
and the inertia tensor of the mass contained in the volume F 

J al = T)oj dqg CT q- 1 , (9) 



■Jr 

are introduced. The dot in equation ([?]) represents the operator d/dt, the symbol d a = d/dq a and e a /3~/ is the completely 
antisymmetric Levi-Civita tensor, £123 = 1. Summation over repeated Greek indices is understood. 

Equation ([?]) shows that the linear angular momentum i/ 1 -* is in general non-zero because the principal axes of the 
inertia tensor Jap, which depend only on the (irregular) shape of the volume V, are not aligned with the principal axes of 
the deformation tensor T>^, which depend on the location of neighbour matter fluctuations (see the discussion in White 
1984). This fact can be illustrated by writing explicitly one component of L in the inertia tensor's principal-axes system, e.g., 
Li oc (^33 — J722): only the off-diagonal elements of the deformation tensor enter. 

In addition, the time dependence of L^ 1 ' is governed by the function a(t) 2 D(t), which simplifies to \t^ 2 t in the Einstein- 
de Sitter universe, as first noted by Doroshkevich (1970) (the convention a = D — (t/to) 2 ^ 3 is adopted). Finally, we emphasize 
the spherical case: if F is a spherical Lagrangian volume, then J a -, = ^ r/oq 5 and 1$ oc e a/ 3 7 2?S = 0; the symbol b~ a j} 
indicates the Kronecker tensor. Consequently, the matter initially contained in a spherical volume does not gain any tidal spin 
during the linear regime (see the discussion in White 1984). Of course, as we will show in the next section, this also holds for 
the variance, i.e., (L^ 2 ere ) = 0. 

Note that the angular momentum may be written directly in terms of the gravitational tidal field (which coincides 
with the shear field in the linear regime) 

4a S "ft - l ( V ' S<1) )<W = (Spdr - i<WV 2 )Vl , (10) 

as 

(t) = a(t) 2 D(t) e aM rf» J al , (11) 
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since the isotropic part of the deformation tensor does not contribute to the antisymmetric tensor product: this justifies the 
use of tidal angular momentum to indicate the angular momentum originated by the coupling between the traceless part of 
the deformation tensor and the inertia tensor. 

Finally, one can express the initial deformation tensor T>^p as an integral in Fourier space: 

-/(^W^' (p)W(pR) . (12) 

Note that, here and hereafter, we are explicitly assuming that the potential field ipi is filtered on scale R by means of the 
window function Wr{(\). The function W(pR) is the Fourier transform of the smoothing function VFjj(q). The filtering of the 
field tpi on an appropriate scale reflects the restriction to the linear evolution of the protogalaxy, since the non-linear coupling 
between different modes is filtered out. 

As far as one considers the evolution of a single collapsing region of assigned enclosing volume F and velocity field, then, 
to calculate its final angular momentum acquired via tidal couplings with the neighbour matter distribution, it is enough 
to supplement the equation (Q) with the time tu at which the perturbation stops expanding and starts to collapse. After 
i.M the angular momentum essentially stops growing, becoming less sensitive to tidal couplings (Peebles 1969). While for a 
general aspherical perturbation the epoch tM is not well defined (see Eisenstein & Loeb 1995 for the case of an ellipsoid), 
if one restricts study to a spherical top-hat region with uniform overdensity 5 = -^(Gnt/tM) 2 ^ 3 , then the time of maximum 
expansion is reached when S = |j(6tt) 2/3 (see Peebles 1980). 

However, one is interested in calculating the tidal galactic spin averaged over an ensemble of realizations of the gravi- 
tational potential random field tpi, rather than in the evolution of the single perturbation. This is particularly important in 
order to compare the results from iV-body simulations against theoretical ones. This programme is carried out in the next 
section. 



2.2 Ensemble average 

We simplify the previous results by considering the expectation value of the square of L, (L 2 )^,. In linear theory this is quite 
simple to do because (L 2 ),/, = (l/ 1 -* 2 ),/,. Here and hereafter, we indicate the linear angular momentum simply by L, and the 
initial gravitational potential by tp = 

From the expression (Q) of L a , we obtain 

(L 2 )^, = a 4 D 2 € a p~f € a /}i 1 i { Ja^ J a t 1 i T>j3 cr T)f)i rT i)^ ! ~ a 4 D 2 e OL fj 1 £ Q /3'y Ja^ Jcr'^i { 'Dpa'Dp/tr')^ ■ (13) 

Note that the correct procedure to compute {L 2 }^ would require complete understanding of how the different realizations 
of ip stochastically influence the boundary of the generic volume F, which in turn determines the inertia tensor J . This is 
a major unsolved problem. To bypass it, we have assumed in the last step of equation (^) that J is uncorrelated with tji, 
thereby simplifying the calculation considerably; the price to pay is a loss of accuracy. This approximation is equivalent to 
preselecting a volume with given inertia tensor and independently assigning the realizations of the field, or, from a different 
point of view, the last step of equation ([H^) would be strictly correct if the random field i/> were to be restricted to those 
realizations of the ensemble which are compatible with the preselected inertia tensor and one would restrict study to those 
objects only. For example, the Hoffman-Ribak algorithm (Hoffman & Ribak 1991; van de Weygaert & Bertschinger 1996; 
Sheth 1995) can be used to set up constrained realizations of the (Gaussian or non-Gaussian) random density field containing 
objects with a given inertia tensor. 

In the context of the Gaussian peak formalism (section 2.3 and Appendix B), it is possible to compute exactly the (small) 
factor by which the neglect of correlation between inertia tensor and the gravitational potential in equation ( |l3| ) leads to an 
overestimate of the rms angular momentum. Nevertheless, we prefer to maintain the present procedure in this section, both 
for its simplicity and because it is applied when calculating perturbative corrections to the linear tidal angular momentum, 
whose distribution deviates non-trivially from the Gaussian one during the mildly non-linear regime (Catelan & Theuns 1996; 
see also the last paragraph in Appendix B). 

The ensemble average in equation (13) can be reduced to an integral over the initial power spectrum P^{k) using 
equation (p^): 

(Pf, a V p , a ,)^= J j0^Pi3P<rP/3iPa> P^{ P )W{pR) 2 , (14) 

where {4>(p)ip(p')}i> = (27t) 3 <5d(p + p')P^(p): the presence of the Dirac function is necessary to guarantee that the density 
power spectrum has the correct units of volume and the conservation of momentum in Fourier space. 

Observing that in the integral (^) the function Pv(p) W(pR) 2 depends only on the modulus p of the wavevector p, and 
that the volume filtered by the window function W(pR) is spherically symmetric, applying the general relation 
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/ dpp a ppPipgF(\p\) = ^r{6 a p6rf + 8aiSps + &aS&pi) / dpp 6 F{p) , (15) 

J sphere J 

which holds for any function F(p) depending only on the radius p, we finally obtain (see Appendix A for a detailed derivation) 
the particularly simple result 

<L 2 )^ = ^a 4 D 2 (fil - 3fi 2 ) a (R) 2 . (16) 

This result holds for any power spectrum -FV(p) and is independent of the statistics of the underlying field. The quantity 
(Jo(i?) 2 is the mass variance on the scale R, i.e. <jq{R) 2 = (2-7T 2 ) -1 J" 00 dpp 6 P 1 p(p)W(pR) 2 . 

The general expression ( |l6| ) is independent of the details of the shape of the boundary surface of the volume V: it depends 
only on the quantities \i\ and jj,2, which are respectively the first and the second invariant of the inertia tensor J a p- The 
latter are defined in terms of the eigenvalues ti , ii and (.3 of the inertia tensor, 

(tlSH+l2+(3, (17) 
(J.2 = Ht 2 + tlt3 + t2t3 • (18) 

(The third invariant is the determinant, 111213.) We stress that the particular combination /1 2 — 3/i2 is null for a sphere, since 
ti = t2 = 13 for a spherical volume T, whereas for a non-spherical volume, fj, 2 — 3/i2 > 0. 

So far we have described the tidal torques acting on a volume centred on a random point in the smoothed density field, 
which does not necessarily enclose a bound protogalaxy. However, galaxies are expected to form around (high) peaks on 
relevant scales in the density field: in the biased galaxy formation scenario, only the density maxima above a given threshold 
(peaks) of the initial Gaussian density field can eventually form galaxies (Kaiser 1984; Politzer & Wise 1984; Peacock & 
Heavens 1985; Bardeen et al. 1986). In the next section, we explicitly assume the volume F to be centred on a peak of the 
underlying smoothed density field. In this case it is possible to take into account exactly the correlation between the inertia 
tensor of the matter contained in F and the potential (and hence also density) field. In the following we will still use the 
notation (L 2 )^, to denote the average over the ensemble of realizations of ip which are compatible with peaks with preselected 
inertia tensor, (L 2 )^, j. The latter corresponds to the average over the (unconstrained) realizations of off-diagonal tidal field 
(see Appendix B). 



2.3 Tidal torques and density peaks 

In this subsection, we specialize the ensemble average (L 2 )^ to the case in which the volume F is centered on a peak of the 
Gaussian density field. Consequently, the combination of eigenvalues fi 2 ~Sfi2 will depend on the parameters characterizing the 
peak's shape. The main difficulty to overcome, however, remains that of estimating which fraction of the matter surrounding 
the primordial peak will eventually collapse onto the peak itself. In other words, the main problem is to estimate the shape 
of the surface boundary of F. This is of course a very difficult and unsolved problem. 

Previous analyses of the topology of the constant-density profiles in the neighborhood of the peaks of the Gaussian field 
8 showed that the isodensity surfaces S c = ^ c o"o, where v c is a critical threshold in units of cro, are simply connected and 
approximately ellipsoidal, at least for sufficiently high density peaks (Doroshkevich 1970; Bardeen et al. 1986; Couchman 
1987). Hence, we assume that the surface boundary of the volume P is determined by the criterion that all the matter above 
the threshold v c — will collapse and form the final galaxy. A similar convention is adopted in Heavens and Peacock (1988). 

A density peak is characterised by the conditions of zero gradient, V5 = 0, and negative definite mass tensor d a dpS. The 
location of the peak is assumed to coincide with the origin q = 0. In the vicinity of the maximum, and in the coordinate 
system oriented along the principal axes of the mass tensor d a dp5, the density field can be approximated by (Bardeen et 
al. 1986) <5(q) = 3(0) — i ~^2 a ^aQa, where \ a , a — 1, 2, 3 are the eigenvalues of the tensor —d a dpS. Since the peak corresponds 
to a maximum, all A Q are positive. Furthermore, it can be assumed that, e.g., Ai > A2 > A3 > 0. Expressing the height of the 
peak in units of <7o, 5(0) = uao, the surface boundary of F may now be written as A^q\ — 1, where the quantities A a 
are the principal semi-axes of the ellipsoidal isodensity surface 5 C = 0: A 2 = 2voo/\ a . Note that this last relation holds only 
if v > 0: therefore, all those objects corresponding to local maxima located in regions of below average mean density (void) 
are not taken into account. Fortunately, such objects will probably fail to form a prominent bound object anyway. 

At this point, noting that the inertia tensor J a p is diagonal in the mass tensor eigenframe, and that F = ^Ai^As, one 
has (see Appendix B for technical details) 

H=\J (L 2 )^ =£C, , (19) 
where we have defined 

£„ = a 2 D r/ a R 5 , , (20) 
and 
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967r :d-7 2 ) 1/2 ^V /2 #C- (2D 



VW \ix J 8(e,p) 3 / 

The meaning of the peak shape parameters, whose distribution V(v,x,e,p) is given in Appendix B, is the following. The 
parameter x is an indicator of the 'sharpness' of the peak, sharper peaks have higher x; the parameters e and p characterize 
the asymmetry of the isodensity profile: e = (Ai — Aa)/2(Ai + A2 + A3) and p = (Ai — 2A2 + Aa)/2(Ai + A2 + A3) . More in 
detail, the parameter e(> 0) measures the ellipticity of the matter distribution in the plane (91,53), while p determines the 
oblateness (0 < p < e) or the prolateness {— e < p < 0) of the triaxial ellipsoid. If e = 0, then p = 0, and the ellipsoid is a 
sphere. The functions A(e,p) and £>(e,p), which are polynomials in e and p, are given explicitly in Appendix B. For a sphere, 
.4(0,0) = and 8(0,0) = 1. 

The meaning of the spectral parameters 7 and R, instead is the following. The parameter 7 is a measure of the width 
of the power spectrum: if the spectrum is a delta function, then 7 = 1; on the contrary, if p 7 P^,(p) is constant over a wide 
range of p, then 7 <C 1; the dependence of 7 on scale for a Cold Dark Matter (CDM) power spectrum is displayed in Fig. 1 
of Bardeen et al. (1986): typically, 7 ~ 0.62 on galactic scale. The parameter _R* is a measure of the coherence scale of the 
density field or, in other words, it gives an indication of the smallest wavelength in the power spectrum: R, is related to the 
smoothing length 7?. For a Gaussian filter and a power-law spectrum, R — yMn+~5) /6 R, , where n is the spectral index. To 
select masses on galactic scale, one typically takes R w 0.5 ft -1 Mpc (for further details see Bardeen et al. 1986; Appendix B 
of this paper). 

The rms angular momentum C in equation jl9| ) has been written as a product of two factors. The first factor £„ = C*(t) 
contains the temporal dependence and can be thought of as the angular momentum unit. On the other hand, the second 
factor £ = £(y, x, e,p; 7) is dimensionless and depends explicitly on the peak shape parameters; in addition, £ depends on the 
power spectrum too, through 7Q The unit C* and the shape term £ are the 'equivalent' of J re /(t) and j e in equations (10) 
and (13) in Heavens & Peacock (1988), respectively. (Specifically, £„ is identical to J re f{t) for / = 0, but £ is averaged over 
the constrained field ip whereas j e corresponds to one single realisation of ij) , £ = {je)l/ 2 -) 

The result in equation ([jl]) is quite simple as a consequence of our factoring out of the dependence of angular momentum 
on the underlying random field ip on the one hand, and on the peak shape parameters u,x,e,p on the other hand. This 
is already manifest in the original equations ( |l3| ) and (^) and allows a simplification of the calculation of the probability 
distribution of £, since V(v,x,e,p) and the second moment of the off-diagonal shear components is all we need to know. 
Heavens & Peacock (1988) calculated the probability distribution function of the modulus L = |L| of the peak's angular 
momentum. Since L depends on the shear components, they needed to compute the joint distribution of shape parameters 
v and \ a (i.e. u, x, e and p), and shear field £ a p; the distribution of early torques on a peak of given height v then leaves a 
four dimensional numerical integration. As we will show in Section 3.2, our averaging over the constrained realisations of the 
gravitational potential narrows only slightly the probability distribution of the angular momentum (see Fig. 7 below). 

In the next section, we will compute several distribution functions for different spin parameters: for this purpose it is 
useful to express the mass M = rjo F enclosed by the isodensity ellipsoid in terms of the peak shape parameters: 

-/ \ 3/2 

M = mhl, = V8 f B(e,p)' 1/2 M* , (22) 

where the reference mass is M* ee rjo i?*. We stress that the mass of the collapsing object is estimated as the mass within 
the volume F: the filtering procedure, which defines the mass unit M», has to be thought of as corresponding to some genuine 
physical filter (more than an artificial one) which prevents the coupling between the non-linear (small-scale) and the linear 
(large-scale) modes in the power spectrum of density fluctuations. The term m(y, x, e,p; 7) describes the deviation (reduction 
or amplification) of the peak mass with respect to the mass unit M* in terms of the shape parameters. We now set the 
smoothing scale R, = 0.23 (0.62/7) 1/2 ft" 1 Mpc in order to have M* = 1.45 (0.62/7) 3/2 10 10 h' 1 M . This ensures that mean 
values of M in equation (j2^) correspond to physical masses ~ 1O 11 M0 typical of spiral galaxies and so allows comparison of 
our angular momentum distributions with galaxy data. To compare against e.g. cluster data, one has to choose an appropriate 
filtering scale i?*. 

It is well known that for individual objects the mass given by ( p^ ) is a rather poor estimate of the mass of the final 
collapsed object, as , e.g., shown by numerical simulations. Indeed, this holds true even for the more sophisticated version of 
the peak formalism like the excursion set method (see the review in White 1994). Even so, the mass- function predicted by 
the latter method, the Press-Schechter function, is remarkably close to the measured mass-function. In the same spirit, we 



* Actually, the factor (1 — 7 2 ) 1 / 2 in equation (21) cannot be obtained directly by specializing equation (16) to the case of Gaussian 
peaks, because it is related to the correlation between inertia tensor J and the random field -0, which has been neglected in equation (13) 
— see Appendix B for a self-consistent derivation. Bearing this in mind, we will continue using the notation {•)ti l , but write the factor 
(1 — 7 2 ) 1 / 2 explicitly in the equations that follow. 
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suggest that the statistical predictions which follow from our adopted procedure are more reliable than the predictions for 
individual objects. 



3 ANGULAR MOMENTUM PROBABILITY DISTRIBUTION FUNCTIONS 



In this section we compute the probability distribution of three spin parameters £p = 1/mP: (i) j3 = 0; (ii) /3 = 1, the specific 
angular momentum, and (iii) /3 = 5/3, which is suggested by the dimensional dependence of angular momentum on mass (see 
below). 

The probability distribution of lp can be computed combining the relation lp = lp(a.i) with the probability distribution 
V(otk), where the at (k = 1, . . . ,4) denote the peak shape parameters. For V(ptk), we will use the distribution function given 
by equation (^) , appropriate for Gaussian peaks. 

A useful statistic is the conditional probability distribution for e.g. lp, given a specific value of one of the parameters a k , 
e.g. given the height of the peak. Such a conditional probability distribution P(£p\ai) dip is formally computed from 



P(£p\ai) dip = Ci 



da 3 / da 4 P(a k ) 



d( n 



(23) 



where Ci is a normalization constant such that J V(£p\cti) dip — 1. The formal expression ( |23[ ) is general, but we will usually 
specialize to the case ai — v and in Appendix C we will show how to compute the distributions for the different cases 
P = 0,1, 5/3. 

The same procedure can be applied to the joint probability distribution of spin £p and mass m, V(m,£p) dm dip: 



V(m,ip) dm dip — / das I doL 4 V{a.k) 



d(m,lp) 



d(ai,a 2 ) 



dm dip 



(24) 



where m(ak) is the mass of the peak in units of M*. Again, we will show in Appendix C how to specialize the parameters a k 
for the different values of f3. 



3.1 Linear growth era 

During the linear growth era, before the peak decouples from the expansion of the Universe, angular momentum i and peak 
mass m are related to the peak properties through equations (|5l]) and (|5^) respectively. Noting that £ and m depend on v 
and x only through the combination y = v/x, we use equation (M) to obtain: 

V(m,£p\v\ < v < v-z) dm dip = 

(/■oo rv-2/y \ 

/ dy W[e(y,m,£p),p(y,m,£p)}J^ p 1 J^ 3 dx x 9 exp[-f(y)x 2 } ) dm d£p . (25) 

Jo Jfi/y / 

This result is obtained by doing the transformation (e,p) — > (m,£p) in two steps. Hence, J ep denotes the modulus of the 
determinant of the transformation (e,p) — -> (A,B) and Jab is the modulus of the determinant of the Jacobian of the trans- 
formation (A,B) — * (m,ip). The function f(y) appearing in the argument of the exponential follows from equation (|o2|). 
The integral over the peak sharpness x can be performed analytically, yet the inversion (A,B) — » (e,p) cannot, as it involves 
finding the real roots of a polynomial of degree six! Further technical details can be found in Appendix C. 

The last equation is a particular case of equation (^) since V{m,£p\v\ < v < V'i) dm d£p is the number of peaks per 
unit volume, in an interval dm dip round specific values m and ip, which in addition have height v between v\ and This 
is because we wish to identify a certain class of objects (e.g., spirals, ellipticals, clusters) assigning a finite //-interval rather 
then with a single value of v (see the discussion in Blumenthal et al. 1984). 

The joint distribution of m and £ (f3 = 0, v\ — 0, vi = oo) is shown in Fig. |l| for 7 = 0.62, typical of the CDM 
model. Note that the most probable value of £, given m, is strongly correlated with m, causing equiprobability contours to be 
stretched approximately along £ oc m 5//3 . This dependence is obvious on dimensional grounds, since, from Eqs. (|2^), £ oc 
and M oc Rt, hence £ oc M 5 ^ 3 . It follows that £/M 5 ^ 3 is independent of the (rather arbitrary) smoothing scale R, (see the 
discussion in Heavens & Peacock 1988). This scaling provides an excellent fit to the actual scaling of I with m and is shown 
in Fig. ^ for comparison. 

The joint distribution of m and the specific angular momentum £/m (f3 = 1) is shown in Fig. ^, again for 7 = 0.62, for 
the whole range of v as well as for restricted ranges, 1/2 < v < 3/2 and 3/2 < v < 5/2. We re-iterate that our description of 
angular momentum for protoobjects is limited to objects with v > 0, hence the total number of objects displayed in Fig. is 
only 96% of the total peak number density for this value of 7 i.e., 4% of peaks have v < for 7 = 0.62. (Note that peaks with 
v < are maxima in a locally underdense region.) The total number of objects with 1/2 < v < 3/2 (3/2 < v < 5/2) in Fig. ^ 
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is 40.1% (37.5%) of the total number of objects shown. The equiprobability contours for peaks with the shown restrictions 
on v display a steeper dependence, i oc m 8 / 3 , however, the most probable £ versus most probable m, averaged over v, again 

5/3 

scales as £ mp oc mm P . The figure also shows that objects with the larger value of v tend to be more massive and in addition 
have a slightly larger £/m as well. 

The dependence on peak height v is further illustrated in Fig. ^| which shows the joint probability distribution V(m,£5 ), 
which, in view of the strong dependence I oc m 5//3 , is the more appropriate quantity to graph. Note that the most probable 
value of l/m 5 / 3 is very insensitive to u, however, the dispersion round the mean is large. 

Conditional probability distributions for £ and m for given peak height v are shown in Figs. ^ and ^| which also illustrate 
the effect of the spectral parameter 7. Higher v peaks have a higher median £. This effect is stronger for lower 7 yet the 
distributions are very wide, as was found previously by Heavens & Peacock (1988) though for the distribution of the modulus 
of the angular momentum. 

Figure ^ shows the distribution of masses, given v. In contrast to Fig. ^j, we have normalised this distribution to 
J V(m\v) dm dv — 0.016/i? 3 , to illustrate the different number of objects of given v, in addition to the distribution as 
a function of mass. Again distributions are wide, with number of objects falling to 1% of the peak value for masses differing 
by factors of five either side of the mean mass. 

So far we have restricted ourselves to the growth of the angular momentum in the linear regime. However, this era of 
growth is not generally prone to comparison with observations (but it can be compared against numerical simulations!) since 
observed structures (galaxies, clusters) have already collapsed. Fortunately, objects presumably acquire most of their angular 
momentum before they collapse, since tidal torques are much less effective afterwards (Peebles 1969, Catelan & Theuns 1996). 
However, the non-linear processes during collapse lead to a redistribution of angular momentum over the different subunits 
that will make up the final object in a complicated way (White 1984; Barnes & Efstathiou 1987; Frenk 1987). 

We will follow Peebles (1969) in assuming that the linear growth effectively ceases after maximum expansion of the object 
and identify the angular momentum at that time with the 'final' angular momentum. This might be a partial description of 
reality. 
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m m 



Figure 2. Same as Fig. |l| but for the specific angular momentum l/m. Left: cquiprobability contours for all objects with v > 0. Right: 
dotted cquiprobability contours correspond to objects with 1/2 <u< 3/2 and continuous contours to objects 3/2 < v < 5/2. Scalings 
(. oc m 5 / 3 and Z oc m 8 / 3 are indicated. 



3.2 Spin distributions at maximum expansion time 

Let us assume that the angular momentum growth described by the equation ([?]) stops as the protogalaxy separates from 
the overall expansion and collapses back onto itself. A fiducial final angular momentum L/ (i.e. £/) can be calculated as the 
value of L at the time tu when, say, 5{M) = —D(tM)V 2 ip = —Dm\/ 2 iP = 1. From equation (Q) and denoting the mass and 
radius of the collapsing object by M and R respectively, we get: 

r 2 A t-,2 , n rr,2 2 Dm T. tr> 2 Dm -2/3 J/f 5/3 Dm (CLm\ 4 ^ 3 ,j5/3 r nc \ 

Lf ~ a M DivfV ipMR = om — — MR oc — — p, J M ' oc — — M 1 . (26) 

L>M L>M L>M \ O-M / 

Since in the Einstein-de Sitter universe a oc D oc t 2//3 , the previous result simplifies to (White 1994) 

L f = L(t M ) <x M 5/3 t 1 ^ 3 . (27) 

This expression will be useful when analysing the dimensionless spin parameter in the next section. Note that, since tM 
depends on M, the latter equation does not imply that t 1 ^ 3 is the full temporal dependence of Lf on tM - see below. (In 
a more general Friedman model, an additional dependence on the density parameter f2 appears, which is approximately a 
scaling Lf oc fl~ ' 07 . Hereafter we will limit ourselves to the EdS case.) 

In the case M refers to the mass of a peak of height v, the maximum expansion time tM depends implicitly on v, which we 
emphasize by writing tu = t M ( v ) = tv- The scaling of t v with v may be recovered extrapolating the linear spherical model to 
the maximum expansion time, giving t v = v~ 3 ^ 2 i„=i . From equation ( pp| ) we find that the 'final' ensemble averaged angular 
momentum contains an extra dependence Cf oc j/~ 3//2 on peak height in comparison to the linear £: 

c s {t v ) = ec»(t v ) = i v - 312 <c»(w) = e f c.(t u =i) . (28) 

The argument goes as follows: start from equation (20) for the unit of spin £*, which we want to compute at the maximum 
expansion time (defined by vD(t v )ao = vD v oq — 1 for a perturbation of height v): 

£ t (U) = alD vVo a Rl = Vo (uD,a )Rl = ^^VoRl- (29) 

vu v vD v 

Consequently, since a 2 , D V ID V oc tl /3 in the EdS universe, 

£,(*„) = v' 1 tl /3 ri R 5 , oc U ■ (30) 
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Figure 3. Joint probability distributions V(m, < v < V2) for CDM 7 = 0.62. Equiprobability contours are drawn at 

log[m £5 V(m, £ 5 ; v\ < u < V2) ln(10.) 2 ] = — 9, . . . , — 3, — 2. Different panels correspond to different ranges [ui,U2], as indicated. The 

33 

dotted contours are identical in all panels and correspond to the whole range [u\ = 0, V2 = °o). 



Evaluating this for v = 1 gives the expression (|2£|)|_[ Explicitly, £„(f„ = i) is calculated at the time of maximum expansion of 
a spherical v — 1 overdensity, D(t„=i)<7o = (3/20)(67r) 2,/3 : 

£.(W) = a^^l^^**!^ 



x 5/3 



45 x 10 10 h^M, 



Here, is the Hubble constant, Ho = 100 /i km/s/Mpc is the present day Hubble constant, z is the redshift. On the 
other hand, £f is given by 

If - -" 3/2 < = -7=f (1 " 7 2 ) 1/2 -" 3/2 (-) 5/2 • (32) 

V15 3 7a; B{e, P y /2 



t The same scaling L* oc 1 M 5 / 3 t^ 3 for the modulus of the angular momentum may be recovered directly from the dimensional 
analysis in equation (Ed), but remembering that (a 7^ /3) Lf ~ dadpip ~ v — 1 V 2 !/ 1 ; we recall that the off-diagonal elements of the shear 
tensor are independent of f while the trace V 2 i/> is proportional to v. 
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Figure 4. Conditional probability for angular momentum £ during the linear regime, given the peak height v, for different values of the 
spectral parameter 7. Different symbols denote different values of v, as indicated in the figure. These plots can be compared against the 
corresponding ones in Fig. f in Heavens & Peacock (f988) - see the discussion in the text. 



The first equality can be compared with equation (16) in Heavens & Peacock (1988), which exhibits the same scaling z/ -3 ' 2 . 

Figure [] illustrates the dependence of if /m 5 ^ 3 at maximum expansion as a function of m for several values of the spectral 
parameter 7. In all cases, equiprobability contours are elongated along if /m 5 ^ 3 oc m _1 , a scaling which follows from the if oc v 
dependence read from equation (|j^) and the dependence m oc v 3 ^ 2 , from equation (|i^). Note that the higher 7 distributions 
are significantly more peaked round their maximum. 

The conditional probability of final angular momentum if for peaks of given height v is plotted in Fig. ^ for several 
values of the spectral parameter 7. Superposed is the corresponding curve (for 7 = 0.5) of Heavens & Peacock (1988). 
The comparison is complicated by the fact that we are plotting in the same diagram the probability distribution for two 
different quantities, essentially V{Lf/C*\v) for the modulus of the angular momentum and V '((L 2 ;}^/ 2 / ' C*\v) for the rms 
angular momentum. This introduces a very small shift between the median values of |L/| and (L 2 )^ 2 which we estimate 
by computing the ratio x = (L/) 1 ^ 2 /(|L/|). In the case of a rotationally symmetric ellipsoid, with Jw = J22, we find 
that x — (D 13 + ^ 23) 1// V((2-'i3 + 23 ta) 1 ^ 2 ) — \/4/7r w 1.13, using the fact that 23 a(3 is Gaussian distributed with variance 
a o (1 — 7 2 )/15- As Fig. (JtJ) testifies, the distribution of £/ is slightly narrower than that of Lf (because of the shear averaging 
already performed) but is still very broad; we stress that the two spin distributions agree extremely well for medium and 
large values of the angular momentum. From this we conclude that the rms angular momentum Cf is a good estimator for 
Lf and that the loss of information, caused by averaging over the constrained realisations of the field, appears to be small. 
These results are useful when discussing the galaxy morphology versus rotational properties. 
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Figure 5. Distribution of mass for peaks of given height, for several values of the spectral parameter 7. Unlike Fig. |4|, mass probabilities 
shown here are not normalised to one, but reflect the fact that there are fewer peaks of higher v. Vertical lines in panels 7 = 0.5 and 
7 = 0.7 denote the values of the most probable mass, given v. obtained from Peacock & Heavens (1990) which fit well. Different symbols 
denote different values of v, as indicated in the figure. 

4 GALAXY MORPHOLOGY AND ANGULAR MOMENTUM 

In this section we compare our predictions against observational data on both the dimensionless spin parameter A and the 
specific angular momentum L/M. 



4.1 Dimensionless spin parameter 

A conventional parametrization of galactic angular momenta is the dimensionless spin parameter A 

\ = L\E\ 1/2 G~ 1/2 M~ 5/2 , (33) 

which is a measure of the ratio between the observed angular velocity u> and the angular velocity u>o which would be 
required to support the system centrifugally (such as e.g. a rotationally supported self-gravitating disc): A w uj/ujq ~ 
[L/(M R 2 )]/[(GM / R 3 ) 1 ^ 2 ], where L is the modulus of the angular momentum, E is the total binding energy of the pro- 
toobject. 

Typically, A depends on the galactic morphological type, being as high as A w 0.5 for spirals and SO galaxies, but only 
A « 0.05 for ellipticals, although the dispersion around these values is large [Efstathiou & Jones 1979; Fall & Efstathiou 1980; 
Kashlinsky 1982; Fall 1983; Davies et al. 1983, Efstathiou & Barnes 1984. More in detail, bright giant ellipticals (magnitude 
Mb < —19) rotate more slowly than faint ellipticals (Mb > —19), and the latter are observed to rotate nearly as rapidly 
as bulges of spirals and lenticulars of the same luminosities (Kormendy & Illingworth 1982; Illingworth & Schechter 1982). 
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Figure 6. Joint probability distribution V{rn,l$ ,), where £5 f = £f/m d ^ 3 , at maximum expansion time for different values of the 

spectral parameter 7. Equiprobability contours are drawn at log[m£s t V(m,£i, ,) ln(10.) 2 ] = — 9, . . . , — 3, — 2. The dashed line in each 

3 3 

panel is the same and indicates the scaling £5 , oc m _1 . 



Also, the flattening of giant ellipticals is thought to be due to velocity anisotropy rather than to coherent rotation (Binney 
1978)]. We will re-examine below the question of the correlation between the galaxy morphology and the amount of angular 
momentum. 

A theoretical estimate of A may be determined from its ensemble average, A = ( A 2 )^, . From equation (^), one 
immediately obtains 

A = C\E\ 1/2 G- 1/2 M~ 5/2 , (34) 

where we have assumed that the total energy E is the same for any realisation of the potential field ip (i.e. it is independent 
of the deformation tensor T> a p). This is of course a questionable assumption. It is well known how difficult in general is the 
exact calculation of E for a generic enclosing volume Y. Hoffman (1988), expanding the total energy to first order in 8, the 
gravitational potential to second order in the position and assuming an ellipsoidal isodensity boundary for F, obtained an 
expression for E in terms of the deformation tensor, E oc 1> a gJ a g. Such a dependence on the deformation tensor of E, when 
consistently taken into account, would hugely complicate the calculation of A. 

However, the order of magnitude of the total energy may be quantified applying the spherical top-hat model to the evolu- 
tion of the overdensity contained in F (see, e.g., Peebles 1980): we expect that such a model becomes a better approximation 
as the height of the peak increases, since high-z/ peaks tend to be more spherical. The spherical model is adopted for the same 
purpose in Heavens & Peacock (1988). 
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Figure 7. Conditional probability distribution for final angular momentum If given the peak height v for several values of the spectral 
parameter 7. Continuous lines: v = 1, dot-short dashed: v = 2, short dashed: v = 3, long dashed: v = 4, for 7 = 0.5, 0.7 and 0.9. 
Lower right panel (7 = 0.5), a comparison against the Heavens & Peacock distribution is shown (see discussion in text). Thick line: 
distribution for u = 1; thin line: fit for v = 1 of the distribution shown in Figure 2a of Heavens & Peacock (1988), V(lj \v = 1) = 
8.2/[(0.89/<?j-) ' 25 + (^^/0.89) 0,47 ] 5 . Note the excellent agreement between these two distributions for medium and large values of the 
angular momentum. 



Since the total binding energy of the protogalaxy scales with maximum expansion time as E oc M 5 ^ 3 t~^^ 3 , then equa- 
tion (^) implies that the final spin parameter A/ (or, equivalently, A/) does not depend on the collapse time tm (actually, 
this should hold provided that the protoobject is isolated from the rest of the universe and that dissipation processes are 
negligible during the linear evolution). 

The resulting expression for the final spin parameter Af depends on the peak shape parameters as follows 



A 



_ 27 w 2 X/2 _! A(e, P y/ 2 

/ 250 1 7 ' V B(e,p) 2 / 3 



(35) 



where, in the framework of the spherical model, the expression for the matter density of the protogalaxy at the maximum 
expansion time t M ^ — t v , p{t v ) = (37r / A) 2 pb{t v ), and the linear prediction of the root mean square density field, D{t u )a^ — 
(3/20)(67r) 2/3 //- 1 , have been used. An analogous expression for A/, in terms of the quantity £5 * discussed in the previous 
section, is 



A/ 



500 



A/2 



3 ,J 



(36) 



which can be compared against equation (28) in Heavens & Peacock (1988). Note that their equation (28) for A/ corresponds 
to a single realisation of if) (but recall that their probability function for A/ does take the shear distribution into account); 
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Figure 8. Left panel: curves denote the conditional probability distributions for A for different peak heights v for a CDM spectrum 
(7 = 0.62). The squares correspond to the fit given in equation (37). Right panel: cumulative distributions for A for the same v and 7. 



however the scaling \f oc v 1 obtained here is the same: indeed, starting from their definition (19) of js , one finds the scaling 
js oc v~' A / 2 from which one recovers the scaling A/ oc in their equation (28). The underlying reason for this is that the 
ensemble average over the shear does not introduce a v dependence. 

Recently, Steinmetz and Bartelmann (1995) attempted to analyze the A parameter in terms of the peak shape parameters 
(note that their definitions of e and p are different from ours). Unfortunately, their statistical analysis is rather crude since 
they do not compute the probability distribution of their e and p and hence are unable to perform the proper average (A 2 ) 1//2 
(compare their equations (44) and (56)). 

The probability distribution function V(Af\u), obtained from equation (pB]), is displayed in Fig.|§] and can be compared 
against Fig. 4 in Heavens & Peacock (1988), bearing in mind that our plots are for the rms spin parameter. Although these 
authors do not show the case 7 = 0.62, 7 = 0.7 is sufficiently close to the CDM value to allow a direct comparison. 

The fit shown in Fig. is 



A s V{A } \v)^T Q 



1 + a 



A, 



Ao(") 



exp 



log : 



2io g e 



(37) 



where T = 0.68, a = 0.02, 9 = 1.12 and A Q (u) = 0.11 v~ x 1 . This fits the actual distribution well for v = 1-4. Note that the 
median Af w Ao(u) has a slightly steeper dependence on v than expected from equation (^). (In Fig. ^ we used the symbol 
A for Af.) 

It can be seen from FigJ^ that higher z^-peaks have on average lower final angular momentum Af. such a statistical anti- 
correlation between A/ (i.e. A/) and the peak height v has been first found by Hoffman (1986; 1988) and Heavens & Peacock 
(1988). We also confirm that the A/-distributions are very broad (factor ~10), and the shift between curves of different A is 
small in comparison with the width of the distribution. Consequently, a given value of A cannot be unambiguously identified 
with a specific value of v. 

We stress the fact that the curve corresponding to V{Af\v = 1) has the median value A/ = J dAj Af V(Af\v — 1) = 0.15. 
This brings the value for the A/ for the halos of the progenitors of spiral galaxies, if assumed to coincide with v = 1 primordial 
peaks on galactic scales (Blumenthal et al. 1984), within a factor of two of the spin parameter of the luminous parts of observed 
spirals. 

It has often been said that energy dissipating cooling processes are necessary to explain the high spin values of spiral 
galaxies (e.g. Efstathiou & Jones 1979; Fall & Efstathiou 1980), and dissipationless iV-body simulations produce objects with 



© 0000 RAS, MNRAS 000, 000-000 



16 P. Catelan and T. Theuns 



A ~ 0.08, too low for spiral galaxies (Barnes & Efstathiou 1987). However, newer TV-body simulations which include dissipation 
(e.g. Navarro, Frenk & White 1995) still fail to produce objects with high A. 

If our method for describing tidal torques is sound, then the importance of cooling processes to boost the value of A from 
the values found from dissipationless TV-body calculations A sa 0.08 to values A w 0.5 found for observed spirals might be 
partially relaxed. Although we understand that cooling must be important to explain why the luminous matter of galaxies is 
concentrated in the middle of the dark halos, it should not be invoked as the only cause of the increment of A for spirals with 
respect to the one for ellipticals, it could be partially explained on statistical grounds. [Furthermore, we estimate (Catelan & 
Theuns 1996) that the inclusion of mildly non-linear leading-order corrections to L increases the value of A over that of the 
purely linear prediction by a factor w 1.3.] 

Our analysis supports the suggestion that, in hierarchical clustering scenarios, the higher v fluctuations will originate 
rather smaller A values, the v ~ 1 fluctuations on galactic scale being statistically associated with spiral galaxies and higher v 
fluctuations with ellipticals (Sandage, Freeman & Stokes 1970; Kashlinsky 1982; Faber 1982; Blumenthal et al. 1984; Hoffman 
1988). 

However it is important to recall that our estimate of A/ is based on an extrapolation of the linear theory and, as discussed 
in Barnes & Efstathiou (1987) and White (1984; 1994), the linear angular momentum is a relatively poor indicator of the final 
angular momentum of a highly non-linear clump of matter (a 1). Non-linear dynamical effects (outward/inward convective 
motion of matter etc.), observed to operate in TV-body simulations, tend to reduce by a factor ~ 3 the final value of L from 
the extrapolation of linear theory [although the total scatter about the mean relation is also a factor of ~ 3: see the discussion 
in Barnes & Efstathiou (1987)]. Clearly, one cannot hope to describe such highly non-linear effects in an exact analytical way. 

4.2 Specific angular momentum 

An obvious advantage of studying the specific angular momentum L/M , rather than the dimensionless parameter A, is that 
the first does not depend on the binding energy E and therefore is not affected by uncertainties associated with dissipative 
processes. In addition, L/M is not plagued by the problem of how to compute (analytically) or measure (observationally or 
from TV-body simulations) the energy of the object. Consequently, L/M is a more robust estimator of the angular momentum 
of an object and is better suited to compare different methods of obtaining L. (In view of the strong dependence L oc M 5 ' 3 , 
the ratio L/M 5 ^ 3 would be an even better discriminant.) 

A comparison of our predictions on L/M against observations is made possible by the data presented in Fig. 1 of Fall 
(1983), who gives values of M (luminous) and L/M for a set of elliptical, Sb and Sc galaxies. This figure shows that the 
visible parts of the sampled spiral galaxies have about 6 times as much angular momentum as elliptical galaxies of the same 
mass. In Fig. |^ we superimpose the data of Fall (1983) on equiprobability contours of V(M,Cf /M) obtained by translating 
to physical units the probability distribution V(m, £i,/), with £\j = if /m, calculated in the previous section. We stress that 
Fall's masses refer to the luminous parts of the galaxies and so should be converted to total masses to allow a comparison. 
Unfortunately, this procedure is not well defined and we have computed the required scaling factor as V2/X, as suggested in 
Fall (1983) for spirals. Taking A = 0.15 for spirals and A = 0.07 for ellipticals, as suggested by our analysis for v m 1 and 
v m 2, we find factors 9.4 and 20, respectively. The smoothing scale is chosen such that the median mass corresponds to the 
typical mass of a spiral galaxy. The agreement between the maximum probability contour and the location of spiral galaxies 
is rather remarkable. However, it is clear from the figure that our prediction of Cf/M does not fully accommodate the data 
of ellipticals which tend to fall outside high probability contours (but see the discussion in the last subsection), yet more 
definite conclusions cannot be drawn because of the evident smallness of the sample. A larger sample of data would enhance 
the fairness of the comparison. The right panel shows equiprobability contours for objects with a restricted range in peak 
height. The location of spiral galaxies is fitted well when selecting peaks with v ~ 1. Selecting peaks with v « 2 does tend 
to shift contours towards the zone occupied by ellipticals, yet falls short of obtaining good agreement. On the whole, it does 
appear that within our description it is not possible to discriminate between the rotational properties of spirals and ellipticals 
solely on the basis of the initial peak height v. 

It is well known that there is a strong correlation between galaxy type and environment, in the sense that ellipticals 
tend to occur in regions of high galaxy density while spirals tend to avoid such regions. Hence, interactions between galactic 
systems, either in the proto or in the collapsed stages, have an influence on the properties of the final observed structures. 
However, even if it were true that ellipticals form from v w 2 peaks, and although it is well known that, within the peaks 
formalism, higher //-peaks are naturally clustered more strongly, the likely direct effects of such enhanced clustering on the 
tidal torques on the peaks themselves are not completely taken into account here. This is best illustrated recalling the fact 
that the off-diagonal shear components are independent of u, or, in other words, the (off-diagonal part of the) shear field 
does not take into account the effect of the stronger clustering of higher peaks. So one is possibly missing one of the most 
basic ingredients which might generate different types of galaxies. Consequently, it is likely that elliptical galaxies are not 
well described by our formalism. We suggest that, with regard to this matter, the model might be improved by extending the 
Taylor expansion - which led to equation (7) - to higher order. 
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Figure 9. Symbols: observational data as taken from Fall (1983); •: ellipticals, □: Sb and o: Sc. Luminous masses from Fall are converted 
to halo masses by multiplying the former by 9.4 (20) for spirals (ellipticals), keeping L/M constant (see text). Left panel: equiprobability 
contours are drawn at log[£ f V(M, C f /M) ln(10.) 2 ] = —5, . . . , — 2. Dashed line indicates the scaling Cf oc M 5 / 3 , which follows from the 

dimensional dependence C* oc MJ . The theoretical equiprobability contours correspond to a single smoothing scale i?*, chosen such 
that M* = 1.74 X 10 11 Mq for h = 0.5. Right panel: same as left panel but for peaks with height 1/2 < v < 3/2 (dotted contours) and 
3/2 < v < 5/2 (continuous contours). 

Our prediction for the dimensional value for Cf/M versus M on the scale of a rich cluster of galaxies is depicted in 
Fig. |lo[ Selecting objects of masses comparable to such structures leads to typical C / /M values one order of magnitude larger 
than for spiral galaxies. Note that the importance of both non-linear interactions and dissipation are less likely to affect this 
prediction than in the comparison with galaxies. 

Our estimate of the specific angular momentum can also be compared against results from numerical simulations. However, 
the status of those simulations seems to be slightly controversial. For instance, Zurek, Quinn & Salmon (1988) find values 
of the specific angular momentum in agreement with spiral galaxies, but too high for ellipticals (e.g. their Fig. 2), yet in 
the calculations of Navarro, Frenk & White (1995, Fig. 10), values of L/M are found consistent with elliptical galaxies but 
too low for spirals. Note that these latter authors have included dissipation and find that a substantial part of the angular 
momentum of the dissipating gas is lost to the surrounding halo. However, even their halos have L/M too low to fall in the 
region of spiral galaxies, making the discrepancy between these sets of simulations even worse. Franx, Illingworth & de Zeeuw 
(1991) have voiced the concern that numerical simulations might lack the dynamic range to be able to describe accurately 
the problem at hand. Having a simulation box of size A and Nyquist wavelength ~ 1/R amounts to changing the effective 
width of the input power spectrum and hence has the tendency to change 7. This effect is illustrated in Fig. |ll|, which shows 
7 for given box size A and smoothing scale R. The bigger box simulation with R = 1 Mpc reaches the asymptotic CDM input 
value 7 = 0.62. Referring to Fig. ^, it is clear that simulations with identical physical spectrum but different dynamic range 
may lead to values of the typical angular momentum differing by factors of a few, especially for rare events. 



4.3 Rotational velocities versus scale 

We can obtain estimates for the rotational velocity from our values of the angular momentum as a function of mass, for 
scales corresponding to galaxies, clusters and superclusters. The angular momentum derived from equation ( [jl] ) gives typical 
values Cf « 1.8 x 10 67 kg m 2 s" 1 for a v = 1 peak, using the median value if = 5.3 appropriate for a CDM 7 = 0.62, for 
h — 0.5, z\ = 3 and M* = 1.84 x 10 1C '/i -1 Mq (corresponding to a Milky Way mass 1.1 x 10 11 Mq), well in line with previous 
estimates for the Galaxy: Peebles (1969) quotes for the Milky Way a value 2.4 x 10 67 kg m 2 s _1 , comparable to the estimate 
of Fall & Efstathiou (1980). This translates to a circular velocity v R = Cf/MR w 140km s" 1 , using M = 1.1 X lO n Af and 
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log(M/M ) 



Figure 10. Equiprobability contours of specific angular momentum versus mass, using a smoothing scale such that M* = 1.16 X 10 Mq 
for h = 0.5 appropriate for rich Abell clusters. Dotted contours refer to all objects whereas continuous contours refer to the subset 
5/2 < v < 9/2. Contour levels are —5, —4, —3, —2.2. These specific angular momenta correspond to rotational velocities 5km s" 1 , on 
cluster scales of 3 Mpc 
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Figure 11. Effective value of the spectral parameter 7 for a CDM simulation of box size A and resolution R. The input CDM spectrum 
is from the fit (7.10) in Efstathiou (1989). 
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R = 20kpc. Applying the same method, we proceed to compute the typical circular velocity of structures on larger scales. 
For a rich cluster of mass ~ 1O 14 M0 and size R ~ 3 Mpc, we find vr ~ 1.5 x 10 4 km s" 1 kpc/3Mpc w 5km s _1 , where we 
used the most likely value of Cj/M read from Fig. ^ for v ~ 3. This very low value is clearly consistent with the lack of 
observational evidence for cluster rotation (Efstathiou & Barnes 1984). For a supercluster of mass 1O 16 M we use the simple 
scaling vr ~ 5 (psc/p c ) 1,/3 (Msc/Mc) 1 ^ 3 w 10 km s" 1 , taking the density of the supercluster to be equal to the background 
density and the density and mass of a typical galaxy cluster as previously. Note that these structures are likely to be still in 
the linear regime, hence this value must be decreased by the fraction t/tja of the age of the Universe in units of the turnaround 
time of the supercluster. Furthermore, the spherical approximation is clearly poor to describe the structures on supercluster 
scales, hence this estimate of vr(sc) must be considered rather uncertain. 



5 SUMMARY AND CONCLUSIONS 

In this paper we re-analyzed the problem of the acquisition of angular momentum by a protoobject, progenitor of a galaxy or 
a cluster, due to tidal interactions with the surrounding matter distribution. This process is of most interest for gravitational 
instability theories of galaxy and cluster formation. In section 2, we reviewed the dynamical description of the spin evolution 
in the version given by White (1984), hence the motion of the matter patch of fluid is followed by applying the Zel'dovich 
approximation. The expression for the linear tidal angular momentum L of a protoobject obtained in this formalism and 
reported in equation (7), contains a combination of the deformation tensor T> a /3 and the inertia tensor J a p of the matter 
inside the Lagrangian volume T, which by definition encloses the collapsing protoobject. Such a combination of the elements 
of the tensors T> and J is zero if Y is a spherical volume (or if the boundary of Y is an equipotential of the gravitational 
potential). The temporal evolution of the spin is governed in the linear regime by the function a(t) 2 D(t), where a(t) is the 
scale factor and D(t) is the growth factor of the linear density perturbations, which is proportional to the cosmic time t in 
the Einstein-de Sitter universe (Doroshkevich 1970). 

These results may be notably simplified by considering the ensemble expectation value of the square of the angular 
momentum (L 2 ). The resulting general (albeit approximate) expression in equation (16) is one of the main results of this 
paper. It highlights the very simple fact that, in linear regime, the ensemble average (L 2 ) is a function of the first and second 
invariant of the inertia tensor J alone and is in addition proportional to the mass variance ao(-R) 2 on the scale R. The fact 
that (L 2 ) depends only on the invariants of the inertia tensor and not on the detailed shape of the surface boundary of Y 
provides a considerable simplification of the calculation. 

In section 2.3, we specialised the statistic (L 2 ) to the case in which the volume Y is centered on a peak of the underlying 
Gaussian density field. The ensemble average is therefore restricted to those realizations of the density field that are compatible 
with the preselected shape of the collapsing object. This constrained ensemble average corresponds to an unconstrained 
ensemble average over the off-diagonal shear (see the extensive discussion in Appendix B). The analysis of the tidal torques 
acting on matter in the vicinity of local density maxima (peaks) during the linear regime has been first attempted by 
Hoffman (1986; 1988) and Heavens & Peacock (1988). In particular, Hoffman first analyzed the correlation between the 
height v of the peak and the amount of angular momentum it acquires by tidal interactions. Heavens & Peacock performed 
a more sophisticated investigation, taking into account also the asymmetric shape of the matter peak, described in terms of 
the eigenvalues — Ai,— A2, — A3 of the mass tensor dadpS (Bardeen et al. 1986). This represents an important improvement 
because the strength of the tidal interactions depends strongly on the shape of the object. Our analysis uses an approach 
complementary to the one of Heavens and Peacock and extends their work. 

The resulting expression for (L 2 ) appropriate for peaks has been recast in equation ( |l9| ) into two factors: £ = (L 2 ) = 
£C*, where I — £(v,x,e,p;y) is dimensionless and contains the dependence on peak shape (height u, sharpness x, ellipticity 
e and prolateness p) and C, = C(t;R r ) is the angular momentum unit which contains the growth rate and is fixed by the 
cut-off R, of the power spectrum. This factorization should be compared with the one in Heavens & Peacock (1988) for the 
modulus J of the angular momentum, their equation (9), which corresponds to a single realisation of the shear field. We 
then proceed by computing the probability distribution of several spin variables (e.g. of i, £/m, £/m 5 ^ 3 versus mass m of 
the peak; m is the dimensionless mass of the peak in units of the mass selected by i?*) by using the probability distribution 
of peak shape parameters (Bardeen et al. 1986). This is an important extension of the work of Heavens & Peacock (1988), 
who analysed the rotational properties of galaxies directly in terms of spin parameters versus the peak height v. The latter, 
however, does not allow to determine uniquely the mass-scale of the progenitor. As we argued at the end of section 2.3, we 
expect our distributions of angular momentum versus mass to be more reliable than the estimates for single objects. 

In sections 3 and 4, we computed and discussed probability distributions for £, l/m and £/m°^ 3 and the dimensionless 
spin parameter A, both in linear regime and at the time of maximum expansion. Our findings can be summarized as follows: 

• linear values of the spin £ have a strong dependence on m: £ oc m 5 ^ 3 , see Fig. [I] 

• for objects of a restricted range in v, this dependence steepens to £ oc m s ^ 3 , see Fig. ^ 

• final values (i.e., at maximum expansion time) £/ of £ scale shallower with m, if oc m 2//3 , see Fig. ^ 
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• median values of If decrease with v for high 7 (i.e. steep spectra) but are rather insensitive to v for low 7 (i.e. broad 
spectra) 

• median values of If decrease with 7 for given v 

• median values for A are \{u = 1) = 0.15 and \(u = 2) = 0.07, comfortably close to the spin of spiral galaxies but 
marginally higher than the value typically quoted for ellipticals 

• dimensional values of Cf/M on galactic scales are in good agreement with measured values for spirals but are too high 
to described ellipticals 

• selecting objects on galactic scales with peak height v ~ 1 gives values of Cf/M typical of spiral galaxies. Selecting 
objects with v ~ 2 moves the position of the most probable Cf/M in the direction of the location of elliptical galaxies in the 
(M, C f /M) plane but the shift falls short by a large margin to give good agreement between values of Cf /M for v = 2 peaks 
and those observed for ellipticals. 

• typical values of the angular momentum of a spiral galaxy are predicted to be Cf « 1.8 x 10 67 kg m 2 s _1 with a 
corresponding circular velocity u CiSp i ra i ~ 140 km s _1 . For a rich cluster, we find a typical value w ClC ius. ~ 5 km s _1 while for a 
supercluster ^c,sup.clus. < 10 km s . 

The reliability of our estimate for the ensemble averaged value A/ = ^/(Ay) depends on the extent to which our approximations 
capture the main features of the dominant processes. We have made the following assumptions: (1) linear perturbation theory, 
extrapolated to the mildly non-linear regime (Zel'dovich approximation), can be used to describe the growth of I until the 
maximum expansion time ; (2) tidal torques spin up the matter until the maximum expansion time of the protoobject and 
are negligible thereafter; (3) the mass of the object can be identified with the mass inside an isodensity ellipsoidal surface 
around the peak; (4) the spherical model can be used to calculate the clump's binding energy E - needed to compute A; (5) 
Gaussian peaks formalism (Bardeen et al. 1986) can be used to compute probability distributions for 1$. 

We briefly comment on the appropriateness of these approximations. The use of the Zel'dovich approximation, very 
powerful in describing the mildly non-linear evolution of matter before shell-crossing, may give only a partial description of 
the evolution of the tidal angular momentum of the clumps themselves. Certainly, it cannot predict the very final stages 
of evolution when clumps merge and interact non-linearly, which leads to the present-day galactic configurations. In fact, as 
shown by the numerical simulations of e.g. Barnes & Efstathiou (1987), Frenk (1987), Zurek, Quinn & Salmon (1988), Navarro 
& Benz (1991) and Navarro & White (1994), the tidally acquired spin may change drastically as objects merge or as angular 
momentum gets transported between core and halo. The merging history of the surrounding protohalo is a key factor in the 
determination of the morphology of a galaxy and the merging processes of dense clumps are associated with substantial loss 
of angular momentum to the halo. Cf - as obtained from the extrapolation of the linear theory - is typically a factor of ~ 3 
larger than the final spin of the non-linear object (Barnes & Efstathiou 1987; Frenk 1987). At present, we see little hope of 
computing theoretically such a reduction factor, which is due to non-linear interactions as well as late infall and dissipation 
processes. Consequently, these processes constrain the applicability of perturbation theory to the period before the maximum 
expansion time (Catelan & Theuns 1996). 

The determination of the surface boundary of the Lagrangian volume T is tricky. This boundary determines the mass of 
the object. We assumed that it can be described by an isodensity contour, which is ellipsoidal when the object is centered on 
a high peak. However, a non-negligible part of the luminous matter might have been captured from the very outer regions 
of the proto-galaxies (Ryden 1988; Quinn & Binney 1992), where the ellipsoidal model is presumably a poor description. In 
addition, the formalism applied here is also based on Taylor expanding the density and the velocity fields around the centre 
of the peak. This expansion breaks down far from the peak's centre and consequently is unable to accommodate late infall 
(see the discussion in Hoffman 1988), a process definitely important during the late stages of galaxy formation (e.g. Navarro 
& White 1994). 

In assuming the spherical model to estimate the clump's binding energy E, we have lost the explicit dependence of E 
on the properties of the underlying linear gravitational potential field ip. Consequently, our calculation could be improved by 
taking this dependence into account (see Hoffman 1988, section II), but this would complicate considerably the computation 
of the ensemble average and is beyond the scope of the present paper. 

Finally, in comparing our predictions against measurements of luminous parts of galaxies, one has to bear in mind that 
dissipative processes, not taken into account here, surely have had a major influence in shaping those objects. 
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APPENDIX A 

In this first appendix, we show explicitly how to derive equation (|l6|). The expression in equation jl4| ) may be concisely 
written as 
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(£> /3ct 2V<t')v> = (<W <W' + <W &<r*' + <W <W) * , (38) 

where we have defined 

a roa 2 

^J^l *>p'P*(p)W<pX) a = £. (39) 
Then, from @ 

(L 2 ) = a 4 _D 2 $ e a /3 7 £ a 0'yl Jai Jo'-f* {Sficr S/3'o' + 5/3/3' 5 arJ i + 5 p a i & (3' a) 

fi D $ (^a/37 ^aP"f f JtT'y' ~t~ £aa~/ ^c(G /r y' Ja ,n /' ~T" ^a/37 ^cto^' Jf3~f') J(T"i 
— a 4 D 2 $ [2 <5 7 y + (5/3<r 5 77 ' — 5/37' ^T CT ) «^/37'] ^7 

= a 4 !) 2 $ [3 J- Q(3 J- Q/3 - ( J aa f] . (40) 

This scalar expression does not depend on the particular coordinate system: let us suppose that the axes coincide with the 
principal axes of the inertia tensor; then 

[Jap) = diag(ti, i 2 , is) , (41) 
thus 3j7a/3 J a /3 — (J aa ) 2 = 2 (n\ — 3/i 2 ), where /ii and /i2 are respectively the first and the second invariants of the inertia 



tensor, namely /ii = Tr(j7) = i\ + t 2 + <>3 an d /i2 = 1112 + ti<<3 + t2<<3- Finally, the equation ( |16| ) is obtained. We stress that 
this result holds for any given Lagrangian volume V whose inertia tensor is characterised by the eigenvalues b a ■ 



APPENDIX B 

In this second appendix, we show more in detail how to specialize the ensemble average (L 2 ),/, to the case in which the volume 
F is centered on a peak of the density field: it turns out that combination of eigenvalues /jf — 3/^2 depends on the parameters 
characterizing the shape of the peak. 

We assume that the isodensity profile 5 C — is the boundary surface of the Lagrangian volume F: such a boundary is 
approximately simply connected and ellipsoidal, at least for sufficiently high density peaks (Doroshkevich 1970; Bardeen et 
al. 1986; Couchman 1987). 

A density peak is characterised by the conditions V<5 = and the tensor d a dp5 has to be negative definite: the location 
of the peak is assumed to coincide with the origin q = 0. In the vicinity of the maximum, and in the coordinate system 
oriented along the principal axes of the tensor d a dp5, the density field is approximately described by the first three terms of 
the Taylor expansion (Bardeen et al. 1986): 

«S(q) = 5(0) - 1 , (42) 

where A Q , a = 1,2,3 are the eigenvalues of the tensor —d a dp5: in such a case, each X a is positive in correspondence of a 
maximum, and it can be assumed that, e.g., Ai > A2 > A3 > 0. Expressing the height of the peak in units of ao, 5(0) = vctq, 
the surface boundary of F may be written as 

+ A + if. = 1 (43) 

A 2 A 2 A 2 ' y ' 

A 2 ^3 

where the quantities A a are the principal semi-axes of the ellipsoidal isodensity surface 5 C = 0: 

Al = ^ . (44) 
At this point, noting that if the tensor d a dp5 is diagonal, the inertia tensor J a p is diagonal as well, 

(Joe?) =diag(n,(, 2 ,(-s) = i^oFdiag^ 2 ,^,^ 2 ) , (45) 
5 



it is straightforward to show that the combination /i x — 3/x 2 reduces to 

.2 o.. _ oil Q 4 c -2 _2 „2 fv\ 5 ( a o\ 5 ^(e.p) 

vS/ \U2/ B(e,p) 3 



M 2 -3 M 2 = 2 11 3 4 5- 2 7 r 2 r ?0 2 (~) (^) , (46) 



where we have introduced the polynomials 

A(e,p) = [p(l+p)] 2 + 3e 2 (l-6p + 2p 2 +3e 2 ) , (47) 
and 

e(e,p) = (l-2p)[(l+p) 2 -9 e 2 ] , (48) 
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and the relation F = ^-AiA2A3 has to be used. 

The parameter x = — cr 2 " 1 V 2 5 = a 2 " 1 (Ai + A 2 + A3) is an indicator of the 'sharpness' of the peak; the quantity cr 2 is related 
to the variance of the field dadpS and it is an element of a set of spectral parameters weighted by powers of the wavevector 
squared, k 2 , 

a h (R) 2 = J^ ^k 2h P^{k)W{kR) 2 . (49) 

It results that (k 2 ) = a\ /a 2 and {k 4 } = <j 2 /a 2 . The parameters e and p characterize the asymmetry of the isodensity profile: 
Ax - A 3 Ar - A 3 (5Q) 



2(Ai+A 2 + A 3 ) 2xa 2 

_ Ai — 2A2 + A3 _ Ai — 2A2 + A3 _ . . 

P = 2(Ai + A 2 + A 3 ) ~ 2xG 2 ' 1 ' 
also 

Ai = ^a 2 x(l+p + 3e) , (52) 

A 2 = iff 2 i(l-2p) , (53) 

A3 = ^a 2 x(l+p-3e) . (54) 

More in detail, the parameter e(> 0) measures the ellipticity of the matter distribution in the plane (qi, 93), while p determines 
the oblateness (0 < p < e) or the prolateness (— e < p < 0) of the triaxial ellipsoid. If e > 0, then spheroids with p = e are 
said to be oblate, while spheroids with p — — e are said to be prolate; if e = 0, then p — 0, the ellipsoid is a sphere: however, 
spherical peaks are strongly improbable (see Bardeen et al. 1986 for more details). 

Finally, the ratio ao /02 in equation (46) may be written in terms of two further relevant spectral parameters, 7 and R t : 

7 = <^> = -4" . (55) 

To (72 

A,eV3-. (56) 
02 

The meaning of these statistics is the following: the parameter 7 is a measure of the width of the power spectrum, since 



((k 2 -(k 2 )) 2 ) _ 1 

— (pp — ( } 

measures the relative spectral width: if the spectrum is a delta function, then 7 = 1; on the contrary, if k 7 P^(k) is constant 
over a wide range of k, then 7 <C 1; the dependence of 7 on the scale for a Cold Dark Matter power spectrum is displayed in 
Fig. 1 of Bardeen et al. (1986): typically, on galactic scale, 7 = 0.62. The parameter R* is a measure of the coherence scale 
in the field or, in other words, it gives an indication of the smallest wavelength in the power spectrum (see the corresponding 
discussion in the main text). 

Inserting the relation ao/<J2 = Rl/Sj in equation (46), and using the general expression (16) for the ensemble average 
(L 2 )^,, and the expression for the volume T, we obtain 



5 

I V " ~ ,,u " u "* \ 7 x ) B{e,pf 



2 12 3?r 2 4-2 2 2 D io/ A(e,p) 



( L )* = ^ ttD »^C| z | wr%5- (58) 



This is an interesting expression: since we first calculated the ensemble average, we have been able to write the expectation 
value (L 2 )^, as a function of the four peak pattern parameters v,x, e,p. Following the suggestion in Heavens & Peacock (1988), 
we can write the previous expression as 



C= ^flj?) =IL, , (59) 
where we have defined 

C, = a 2 D 770 a Rl , (60) 
and 

96^ ( v\ 5/2 A{e,pf/ 2 



V15 3 \"ix J B(e,p)' A / 2 

The significance of the quantities I and C, is discussed in the main text. 

The probability distribution V(u,x,e,p) du dx de dp of the peak shape parameters has been obtained by Bardeen et 
al. (1986): it reads 



© 0000 RAS, MNRAS 000, 000-000 



24 P. Catelan and T. Theuns 



V(v,x,e,p) — Vox W(e,p) exp 



2 2 X [6e +P ' 2(1 - 7 2) 



(62) 



in the allowed domain (a triangle) < e < \ , —e<p<e and \ < e < | , 3e— 1 < p < e , and by definition vanishing outside 
of the triangle. The normalization constant is Vo = 3 2 5 5/2 /(27r) 3 i?» (1 — 7 2 ) 1/2 and the function W(e,p) = e(e 2 — p 2 )B(e,p). 

We stress the fact that the multiplicative factor (1 — 7 2 ) 1 / 2 appearing in equation (21), but still missing in equation (|6l|), 
can not be obtained by the procedure outlined so far. The origin of this factor is as follows. The procedure outlined so far 
is equivalent to preselecting a volume V with a given inertia tensor and independently assigning the realisation of tp, hence 
neglecting the fact that the inertia tensor is in general statistically correlated with ip. To obtain consistent results when 
studying objects with a given inertia tensor, the ensemble average of L 2 should not be calculated over the realizations of the 
unconstrained potential ip (as we have done so far - see section 2.2), but only over those realizations of ip which are compatible 
with the preselected F. In general, this is a hugely difficult calculation to perform. Fortunately, in the framework of the 
Gaussian peak model, this can be done both consistently and easily, in that the angular momentum is determined by the 
off-diagonal components of the shear field T> a p only, and the latter turn out to be statistically independent of the shape of the 
peak (i.e., of its inertia tensor). Consequently, as far as the angular momentum is concerned, to compute the ensemble average 
over the constrained potential field, it is enough to compute the ensemble average over the off-diagonal (unconstrained!) 

shear components, {■)t> &- 

More in detail, the joint distribution of the shape parameters of Gaussian peaks and the off-diagonal shear components 
is given by (see equation (A24) in Heavens & Peacock 1988): 

V{v,x,e,p,Vi2,V w ,V 2 z) = T{v,x,e,p) T(V 12 ,V i:i ,V23) , (63) 
where the T> are Gaussian distributed with zero mean and variance 0% = cr% (1 — 7 2 )/15: 

V(V 12 ,V 1S ,V2 S ) = (2ttoI)' 3/2 exp[-(D 2 2 +0 2 3 +2?l3)/2a|,] . (64) 
S8J) we have computed the ensemble average over the (unconstrained) potential field ip, leading to a 



Note that in equation ( 
dispersion 



(2? /3 CT £> p> a >)^, = {Sp a Spiel + Sppi & aa i + Sp a i 5 pi CT ) 



15 



(65) 



whereas computing the ensemble average over the off-diagonal shear (which is equivalent to computing the average over the 
constrained potential field) we now find: 

(l-7Vo 2 



(Dp a T)pi a i)^\j = {T)p tT Vpi tJ i)'D od = {Sp a Spi a i + Sppi S atJ i + 5p a i 5pi a ) 



15 



(66) 



Repeating the calculation in Appendix A, but using $ = a\~,, one recovers equation (^lj) instead of equation (|6l|). 

Intuitively, the smaller dispersion in the latter case makes sense: the random field is constrained to have a peak with 
specified shape, hence the dispersion of the shear tensor will be smaller than in the first, unconstrained case. This is the 
origin of the extra factor (1 — 7 2 ) 1//2 appearing in equation (^): it properly takes into account the reduction in rms angular 
momentum because the reduction of the dispersion in shear due to the imposed constrained of having a peak with given 
shape. 

A final remark: although the procedure outlined in section 2.2 is only approximate, it does have some important advantages 
when one tries to address fundamental issues related to the tidal torquing of the progenitors of the present structures. For 
instance, (i) how does the galactic tidal spin evolve during the non-linear regime? It is know, in fact, that even if the 
primordial fluctuations are Gaussian, the non-linear evolution will ensure that the mass-density fluctuations become highly 
non-Gaussian. The simple Gaussian peak model, for example, cannot be applied during the non-linear regime, yet the present 
approximate method can be applied (Catelan & Theuns 1996); (ii) how does the galactic tidal spin evolve during the non- 
linear regime starting from intrinsically non-Gaussian initial conditions? Again, the simple Gaussian peak model cannot be 
used (work in progress). The procedure suggested in section 2.2 could give some insight into these very complicated issues. 



APPENDIX C 

In this Appendix we give some technical details on how to compute some probability distribution functions, e.g. the one in 
equation (|25| ). We start from the probability distribution V(f,x,e,p), given in equation (162]). The distribution P(m,£p; v\ < 
v < V2) can be obtained from that by transforming (e,p) — > (m,£p). Technically, this can be done more easily in two steps, 
(e,p) — ► (A,B) and [A, B) — * (m,£p), by using the quantities A and B as intermediate variables. To do this we first invert 
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the relations £p = £p(A, B) and m — m(A, B), using equations (fell) and u22|), to find: 



A = JfeJ w (67) 



B - (-)"(-) 3 ' (68) 

where and mo are numerical constants. One proceeds by computing the determinant of the Jacobian of the transformation, 
(A,B) — > (rn,£p), which is Jab = m£p/4AB. On the other hand, J ep = 216e(e 2 — p 2 )(l — p 2 — 3e 3 ) is the determinant of the 
Jacobian of the transformation (e,p) —> (A, B). The parameters e and p occurring in J ep have to be written in terms of A and 
B, but this inversion (e,p) — > (A,B) cannot be performed analytically, as it involves finding the real roots of a polynomial of 
degree six! Combining these results, we find for the probability distribution: 

V(m, lp\ v\ < v < 1/2) dm dip — [ du dxV[v,x,e(is,x,m,£p),p(i/,x,m,£p)] J~ p J^ g dm dip . (69) 

One now observes that the inversion (A, B) — > (m, £p) only depends on the ratio y = v jx of the two variables v and x. Hence, 
the foregoing expression can be simplified further in terms of y to give equation (25) quoted previously. The function f(y) 
appearing in the argument of the exponential in that equation follows from equation (pi^), f(y) = y 2 /2 + 5(3e 2 +p 2 )/2 + (1 — 
7y) 2 /2(l — 7 2 ). One can now perform the integral over the peak sharpness x analytically, leaving a one dimensional numerical 
integration to obtain V(m,lp\v\ < v < 1/2)- 

In some cases (e.g. to compute A or If) there is an extra v dependence which makes that A and B do not depend on the 
ratio u/x alone. In the computation of the distribution V(m,lf /m 13 ), we eliminated v and x in favour of m and If/m 13 , and 
V follows from a two dimensional integration over e and p. The computation of the A distribution for given u, on the other 
hand, can again be obtained from a one dimensional numerical integration, by transforming (e,p) — > (m, A) and noting that 
the latter change of variables only depends on y = mx z ^ 2 and x, which allows one to compute the x integral analytically. We 
do not give the latter analytical expressions since they are rather cumbersome. 
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